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A method is presented that allows exact calculations of fragment multiplicity distributions for a canonical ensemble of non- 
. interacting clusters. Fragmentation properties are shown to depend on only a few parameters. Fragments are shown to be 
^\ ' copiously produced above the transition temperature. At this transition temperature, the calculated multiplicity distributions 
broaden and become strongly super-Poissonian. This behavior is compared to predictions from a percolation model. A 
corresponding microcanonical formalism is also presented. 
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I. INTRODUCTION 



> 

\^ ' Heavy ion collisions where excitation energies are of the order of 10 MeV per nucleon probe the energy regime where 
the nuclear liquid gas transition is expected to take place. Below energies of approximately 50A MeV, symmetric 
collisions are expected to produce sources that evaporate particles as would be expected from a hot liquid drop, whereas 
above this threshold, the excited source is expected to explode, producing larger clusters through simultaneous multi- 
fragmentation. In this energy regime, the process of fragment production is not clear, and comparisons with data 
have been made with a disparate set of models, ranging in simplicity from percolation descriptions |jL^| and lattice 
gas models [|| , to evaporative models pj , dynamical simulations |||]j7|,||,[l(| , and microcanonical samplings Jll],[l^] • 

Fluctuations behave in a special manner in the region of phase transition, so it should seem that the study of 
fluctuations of fragmentation observables might prove insightful for investigating multifragmentation. Moretto and 
collaborators |p| have measured multiplicity distributions of intermediate-mass fragments (IMFs) as a function of 
excitation energy for a variety of projectile/target combinations utilizing beams with energies up to 60 A MeV. The 
Ch ■ analysis showed sub-Poissonian multiplicity distributions that could accurately be described with binomial distribu- 
tions. These observations have inspired a variety of explanations p^-p7|. 

Recently, Chase and Mekjian have discovered a method for exact calculation of the canonical partition function 
for non-interacting clusters. In this paper we extend this approach to include fragmentation observables. We present a 
method for exact determination of both multiplicity distributions and their moments. When raising the temperature 
while keeping the volume fixed, we observe a sharp transition for fragmentation at the same point where the specific 
heat peaks. At this threshold, the multiplicity distribution becomes remarkably wide. We associate this behavior with 
a first order phase transition and remark that the broadening might well disappear in a microcanonical treatment. 

We parameterize the width of the fragment multiplicity distribution, relative to the mean, with a correlation 
coefficient £ which is described in the next section. The method for calculating fragmentation observables from a 
canonical ensemble is presented in section III while the results are presented in section IV . In section [y] we contrast the 



results of this model with those of a percolation model. Expressions that could be used for microcanonical calculations 



are presented in section VI 



II. MULTIPLICITY DISTRIBUTIONS AND THE CORRELATION COEFFICIENT 

A super or sub-Poissonian multiplicity distribution is one whose variance exceeds or falls below it's mean respectively. 
The difference of the variance and the mean can also be written as a correlation. We demonstrate this by considering 
the emission into an arbitrarily large number of states i, each of which is infmitcsimally probable. A state could be 
defined as a specific type of IMF emitted into an arbitrarily small bin in momentum space. The difference of the 
variance and the mean is 
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If the bins are arbitrarily small, one may discard the terms in the second sum proportional to < Hi > 2 , then use the 
fact that nf — rn for n^O or 1, to discard the remainder of the second term and obtain, 

= a 2 - (n) _ Tn^jjiPi - ("»))K- - ( n j)) , - 

4 - (n) 2 ~ (n) 2 1 j 

Here, £ can be interpreted as a correlation coefficient. It is positive if the emission into two different bins is positively 
correlated. The only assumption going into this derivation is that the bins may be divided into arbitrarily small sizes. 

The simplest examples of sub- and super-Poissonian multiplicity distributions are the binomial and negative binomial 
distributions. The binomial distribution is defined in terms of two parameters, p and N. 

N* 

Pn= UN ' ,. p n (l-p) N -\ (3) 

n\(N - n)\ 

where p is the probability of success in one of N tries. For the binomial distribution the mean and correlation 
coefficients become 

<n)bin. = pN (4) 
£bin. = -1/N, (5) 

and stays negative. Like most correlations it is proportion to the inverse of the system size. 
The negative binomial distribution is also defined by two parameters p and N. 

= (N + n-l)\ r ( ,, 
(N -l)\n\ (l+p) N + n ' K ' 

The correlation coefficient in this case is opposite to that of the binomial distribution. 

(n)neg.bin. = pN (7) 
Cnog.bin. = 1/N. (8) 

Binomial and negative binomial distributions result when one considers populating N quantum levels with fermions 
or bosons respectively with p representing the average population of each level. 

Random emission from a large number of uncorrelated sources leads to a Poissonian distribution. The binomial 
distribution suggests that conservation of particle number would give a negative correlation coefficient of order 1/N, 
where N = A/ a is the number of IMFs of characteristic size a that could fit into the system. For sufficiently small 
systems, this negative contribution from particle-number conservation dominates, with the extreme case being where 
a is more than half the system size meaning that no more than one IMF can be emitted. Other negative correlations 
are expected due to energy conservation. If IMF emission requires energy, e.g. escaping a Coulomb barrier, energy 
conservation is expected to reduce the probability of emitting a second IMF. We will see in the next sections that 
positive correlations can arise, principally due to surface considerations. 

Aside from the size of the correlations, it is also of interest to understand whether the entire multiplicity distribution, 
is well described by a two-parameter fit to a binomial or negative binomial distribution. If the reduction of the results 
to two parameters is valid, comparison of different models and data is greatly simplified. 



III. CANONICAL DISTRIBUTIONS OF NON-INTERACTING CLUSTERS 



Chase and Mekjian have shown that canonical partition functions can be easily calculated in terms of the partition 
functions of single clusters. This allows the calculation of thermodynamic quantities for a system of fixed nucleon 
number without resorting to numerically intensive Monte-Carlo procedures. If the partition function for a single 
cluster of size au is denoted by uik, the partition function for a system of size A may be written, 



n k 

(Sn fc a fc =A) fe 

Thus, S7^4 is expressed recursively in terms of u>k and Qyj/ for A' < A. Proof of this relation is given in the appendix. 
The only shortcoming of this approach is that explicit (not mean field) interactions between fragments are ignored. 

Moments of the multiplicity distribution may be expressed in terms of the partition functions. The moments can 
then be used to derive the correlation coefficient, £ a defined in Eq. (Q) or the multiplicity distribution as discussed 
below. The first moment is the mean which is defined as: 

{nh)=Uh ^r (10) 

Rather than considering moments, (nV 1 ), it is more convenient to consider factorial moments, Fb,A,m> defined as: 

Fb,A,m = (nbim - 1) • • • (n b - m + 1)) (11) 
Calculation of the factorial moments for n& defined in Eq. ( [TT| ) is simple if b refers to single species k: 

fk,A,m — W fc . (12) 

' 'a 

However, if b refers to a set of species, rib — Tikeb^ki where the various species that comprise b have different masses, 
Eq. ( p^ ) is no longer valid. One must then generate the factorial moments using the recursion relation, 

Fb,A,m = ^^UkFb,A-ak,m-l ^ °\ (13) 

keb A 

which is true in general. The proofs of Eq. (|l^) and Eq. ( |l3| ) are given in the second section of the appendix. 

As shown in the third subsection of the appendix, factorial moments are sufficient for calculating the entire multi- 
plicity distribution via the relation, 

P b ,A,n = £ F b , A<m - !_(_!)— « (14) 

' — ' (m — n)\n\ 

m>n 

where Pb A,n gives the probability of viewing n fragments of type b in a system of size A. 

More directly, one may also generate the multiplicity distribution, without knowing the factorial moments, through 
the recursion relation, 



(15) 
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Proof of Eq. ( |l5| ) is presented in the appendix. This direct method of producing the multiplicity distribution has 
proven more numerically stable than generating the distribution from the factorial moments. This improvement can 



be traced to the alternating signs in Eq. (14) 



Summarizing the technique, one starts by calculating partition functions for individual fragments ujk- One may 
then generate partition functions, 57^, by using the recursion relation, Eq. (^). The recursion relation for factorial 
moments, Eq. ([l3|), then allows one to generate the factorial moments, which in turn allow the determination of the 
entire multiplicity distribution using Eq. (^) . Alternatively, one may calculate the multiplicity distributions directly 
using Eq. (|15|). The obtained multiplicity distributions are exact. Although the sums used in the recursion relation 
are performed numerically, they require only a fraction of a second of computer time. 



IV. A LIQUID-DROP PICTURE FOR INDIVIDUAL CLUSTERS 



For our purposes, we consider the partition function for individual fragments of mass k as, 



„ 4 = y{2^} 3/ V„.,/r, (16) 

Fk,int = A«* + /.of + /ojt 5 ' 3 , (17) 

where the volume of the system is V, the mass of a single nucleon is M , and the fragment's internal free energy Fk,int 
is split into a bulk term and a surface term. One sees that the bulk term is irrelevant in determining fragmentation 
observables since it factors out of the partition function. Thus, aside from the system size A, all fragmentation 
observables are determined by three parameters, the ratio of the surface term to the temperature f s /T, the intrinsic 
entropy, s = (V/A)(MT) 3 ^ 2 , and Coulomb term, f c . This implies that many details of a system's microscopic makeup, 
e.g. Fermi vs. Bose nature of the internal excitation, are irrelevant in determining the statistics of fragmentation. If 
the surface term is negligible fragmentation is determined purely by the intrinsic entropy. 
For the surface and Coulomb terms, we use the parameters of the nuclear liquid-drop model, 

f s = 17.2MeV (18) 
fc = 0.70 ^1 - (j-^j 1 \. MeV, (19) 

where the form of the Coulomb term was taken to account for the screening of the Coulomb repulsion by the nuclear 
medium in a Wigner-Seitz like parameterization [fl2[ with po referring to nuclear saturation density of 0.15 frn -3 . 

All the calculations presented in this paper assumed a density of one third of nuclear matter density. The behavior 
at different densities is not qualitatively different, with the exception of the relative importance of the Coulomb term. 
The fragmentation transition described below occurs when s — (l/p)(mT) 3 / 2 is of order unity. Therefore a change of 
density affects the temperature where fragmentation sets in. An excluded volume could easily be incorporated into the 
problem by replacing V with V(l — p/po)- This is equivalent to changing the density, and does not qualitatively affect 
the results. The surface energy is chosen to be a constant, / s = 17.2 MeV. One could imagine scaling f s as a function 
of p or temperature, although one might object to incorporating a temperature dependence that is not of the nature, 
e E / T . The choice of f a does affect the transition temperature and it's width. Larger choices of f s lead to sharper 
transitions. The surface term should disappear at the critical point. Although it is easy to incorporate excluded 
volume effects and a running surface term, we choose not to, as our principal goal in this study is to understand the 
general properties of this approach. 

In Figure 0, the mass distribution dN /da is shown for three temperatures, T=7.0, 7.8 and 8.6 MeV. The overall size 
of the system was chosen to be 250 and Coulomb effects were neglected in these calculations. The mass distribution 
has been multiplied by a to emphasize how the composite nucleons are partitioned into the various sized drops. At 
7 MeV, the nucleons are largely contained in one drop, while at the higher temperature, nearly all the particles are 
part of small clusters. 

The average multiplicity of intermediate-mass fragments (5 < a < 40), denoted as IMFs, are shown in the upper 
panel of Figure | for the cases where the density is 0.05 fm 3 and the overall system size is again A = 100. Results are 
plotted against the temperature, both for the case where the Coulomb term is included as well as for the case where 
it is neglected. The inclusion of Coulomb pushes fragmentation down towards lower temperatures. The trend would 
strengthen if we were to consider larger systems. When the excitation energy exceeds the fragmentation threshold, 
average IMF multiplicities quickly climb to over a half dozen per event. 

Correlation coefficients are shown in the lower panel of Figure || with the system size is set at 100. Coefficients are 
shown both for the case where Coulomb is included and for the case where Coulomb is neglected. 

When Coulomb is neglected super-Poissonian behavior ensues at lower excitations as evidenced by the positive 
values of £. At low temperatures, this behavior may be understood by considering the surface penalty for emitting a 
single fragment of mass a, 

T(A^ A -a) ocexp j-a A(i? surf /T)| (20) 

If a second IMF is emitted, it feels the same penalty, except for the fact that dF/dA is evaluated at a smaller overall 
size, A — a. The correlation thus becomes, 

_ T(A — > A — a)T(A -^A~a) 
W ~ T(A - a -» A - 2a) 2 1 ' 

^ a 2 rf 2 -F sur f 
~ T dA 2 
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The dashed line represents the surface contribution to £ and closely follows the statistical calculation at lower tem- 
peratures. The average size of an IMF was used for a. 

The inclusion of Coulomb reduced £ to negative values as seen in Figure 0. The Coulomb contribution to £ may be 
approximated in a similar fashion as the surface contribution, 

U, « ~%^A-V\ (22) 

Again, the approximation is plotted as a dashed line in the lower panel of Figure |[ 

The simple estimate of £ works well at low temperature, where most of the particles reside in a single fragment. 
The approximations for £ improve for larger systems when Coulomb is ignored but become worse for large systems 
when Coulomb is included. This arises because a large system does not wish to form a single drop when Coulomb is 
included. 

At high temperatures the sub-Poissonian behavior may be roughly understood as arising from particle conservation. 
As seen in section |l| when considering binomial distributions, one expects a negative correlation of order l/N, where 
N is the maximum number of IMFs that fit into the system. 

The super-Poissonian (or nearly super-Poissonian when Coulomb is included) behavior at the fragmentation thresh- 
old is especially interesting. Perhaps, this may be understood by stating that the partition function is sampling two 
competing configurations, one with one large fragment surrounded by gas and a second one where numerous IMFs 
are present. At the fragmentation threshold, where the system is undergoing a transition, both configurations occur 
leading to a broadened multiplicity distribution. Figure [| displays £, again scaled by A, as a function of temperatures 
between 7 and 9 MeV, for 3 sizes, A = 250, A — 500 and A = 1000. Coulomb is ignored in these calculations. As 
the system size increases a singularity in £ develops at the boiling temperature. This is related to a discontinuity 



in the energy vs. temperature at constant volume |19|, which is characteristic of this model but not characteristic 
of a typical liquid-gas transition. In a liquid-gas transition C p becomes singular but not Cy. Since the peak in £ 
vs. T is linked to the discontinuity in the energy density at the same temperature, the peak might disappear in a 
microcanonical treatment. 

The multiplicity distributions for the A = 250 system with Coulomb ignored are shown in Figure |] for three 
temperatures, 7, 8 and 9 MeV. They are compared to negative binomial and binomial distributions respectively 
where the two parameters are chosen to fit the mean and variance of the distributions. At 7 and 9 MeV, the two- 
parameter fits seem quite sufficient to describe the multiplicity distributions, whereas at 8 MeV, where the distribution 
is strongly super-Poissonian, the distribution's shape is poorly described by a negative binomial fit. This emphasizes 
that two classes of events, corresponding to the two phases, contribute at this temperature. 



V. COMPARISON WITH THE PERCOLATION MODEL 



Here, we present a brief comparison of the results of the canonical model described in the previous sections with 
a percolation calculation ||. For the percolation calculation, a cubic lattice of sites within a specified radius of the 
center are attached by bonds. The bonds are then broken with a probability p. When p is in the neighborhood of 
p c = 0.7512, the system undergoes a second order phase transition and the size of the average largest cluster rapidly 
falls. 

In the process of comparing percolation calculations to nuclear data, Li et al. p0| have developed a simple mapping 
of the bond breaking probability with the temperature. By using this mapping we are able to plot the results against 
temperature instead of p, and more readily compare to the results of the statistical model presented earlier. In the 
calculation presented here, the number sites in the spherical lattice was 203 and an IMF is defined was a fragment 
with mass between 6 and 40. 

The mean IMF multiplicity and correlation coefficients are displayed in Figure ||. They are contrasted to the 
results of the canonical-ensemble calculations where A was chosen as 200 and Coulomb was ignored. There are three 
striking differences in the results. First, the percolation calculation does not yield as many fragments as the statistical 
calculation, which produces over twice as many fragments. Secondly, the percolation calculation does not exhibit 
a spike in the correlation coefficient as plotted against the temperature. This is expected as the percolation model 
does not contain a first order phase transition. Finally, the correlation coefficient is much larger for the case of the 
percolation calculation at low temperature. This is due to the fragmentative nature of percolation, where the existence 
of a fragment opens surface for the production of a second fragment M . 



VI. MICROCANONICAL CALCULATIONS 



A microcanonical approach, which considers configurations only at a specific energy, would be more realistic since 
nuclear collisions do not take place with contact to a heat bath. The importance of performing microcanonical 
calculations is emphasized by the existence of the first order phase transition, which collapses a wide range of energy 
densities to a narrow range of temperatures. In this section we present expressions for calculating fragmentation 
observables within a microcanonical context. 

It is straightforward to obtain the needed microcanonical quantities from the expressions for partition functions 
by Fourier transforming corresponding canonical objects over a range of complex j3 = 1/T. For instance, the level 
density, p(E), may be obtained from the partition functions through, 

p(E) = Tr a S(E-E a ) (23) 

= -L J df3VL A {i[3)e^ E . 

One may calculate recursion relations for the factorial moments and for the multiplicity distributions at fixed energy. 
First we consider, fb,A,m{i(3) and Pb,A,n(i[3) which are the canonical quantities without the normalization brought on 
by dividing by Q. As described in the appendix, they may be found via recursion relations. 

f b ,A, m (P) = n A (P)F b ^ m (P) (24) 

= ^ Ukfb,A-q h ,m-l{0) 
keb 

Pb,A,n{P) = Sl A {f3) PbA , n {(3) (25) 

= —/UkPb,A-a k ,n-l{P) 
n f— ( 
keb 

The factorial moments and multiplicity distributions at fixed energy E may then be expressed as, 

Fb ^™w - — — ( } 

S df} Pb , A , n {ifi)e ifiE 
Pb ^ n{E) ~ 2MB) 

However, the integrations over (3 needed to obtain the relevant microcanonical quantities do make the calculations 
significantly more numerically intensive. It is not clear to what degree the stationary-phase approximation |2l| might 
allow one to avert the numerically costly integration, especially given the first order phase transition which causes 
discontinuities as a function of [3. 

If one is interested in the microcanonical distribution, for a range of energies in the neighborhood of E, rather than 
for an exact value of E, one may replace the delta function in Eq. (^3|) by a Gaussian, 

5(E-E a )^^L=e Xp ~ {E ~ E 2 a) \ (27) 

V27T77 2 V 

where r\ defines the width of the neighborhood. One may incorporate the broadening of the S function by modifying 
the phase factors used in the Fourier transforms, 

e^-^exp-irfiS 2 (28) 

Numerical implementation of the Fourier transform simplifies for broader widths 77, as it effectively narrows the 
required integration range for (3. 

An alternative way to approach the constraint of energy conservation is to discreetize the energy and measure it 
with integral values. One can then treat energy in the same manner with which one would treat other conserved 
charges. For instance, energy might by measured in steps of 0.1 MeV, with an integer Q measuring the energy. If 
u>i counts the number of states with energy qi and mass <Zj , the number of states of the system of mass A with net 
energy Q becomes, 



N(A,Q)= J2 n ^T 



(29) 



One may then derive the recursion relatio n in a manner similar to the derivation of the recursion relation for the 
partition function shown in Section VIII A . Since there are two conserved "charges" , two recursion relations may be 
derived, 



N(A, Q) = ^2 ^UiN(A -a u Q- ?i ) 



(30) 

(31) 
(32) 



In a manner similar to the derivations in Sections VIII B and VIII D one may derive recursion relations for the 
factorial moments and multiplicity distributions, 



,m — — ak.Z — Qk,m—l 

k£b 

Pb.A.Q.m — y ^fe — Pb. 



;A-A k ,C 



k,n-l~ 



N(A-a k ,Q-q k ) 
N(A,Q) 

N(A-a k ,Q-q k ) 



keb 



N(A,Q) 



(33) 
(34) 



There are two practical differences between these expressions and the corresponding canonical expressions. First, 
an extra index has been added that denotes the energy of the system. Secondly, the index k refers to a set of states 
specified both by mass and charge. In practice this leads to a longer calculation by a factor of the number of energy 
steps squared. Thus if one wishes to perform a microcanonical calculation with an excitation energy of one GeV, 
using energy steps of 1.0 MeV, the required computer time would be expected to increase by 10 6 as compared to a 
microcanonical calculation at a single temperature. 



VII. CONCLUSIONS 



The relations derived in this paper allow the exact calculation of all fragmentation observables within the context of 
a canonical ensemble of fragments. The relations allow the exact summation over all possible partitions of A nucleons 
within a fraction of a second of CPU time. The drawback of this approach is that interactions between fragments 
may only be included in a mean-field context, or through excluded volume. The advantage of this type of approach 
goes beyond numerical convenience. Unlike many Monte Carlo approaches, this technique has no arbitrary choices 
inherent to the algorithm. Thus, when one has finished with a calculation, one can more clearly state what expressions 
one has solved. 

This approach can be made more realistic in three ways. First, the sum over intrinsic states, which here was 
performed in a liquid-drop context, can be replaced by a sum over true nuclear states without a significant cost in the 
time of the calculation. Secondly, isospin conservation may be included. It is straightforward to extend the recursion 
relations to two conserved charges p2j without significantly increasing the complexity of the approach. Finally, 
a microcanonical ensemble would represent a significantly more realistic description. Given that an entire region 
of energy densities is described by a single temperature, the spikes observed in our calculations of the correlation 
coefficient, and in previous calculations of Cy, are expected to disappear in a microcanonical approach. As 
discussed in the previous section, the methods presented here might be extended to a microcanonical description, 
although computational pitfalls may complicate numerical implementation of the expressions derived here. 

We finish by stating our belief that, although this type of approach is not yet to the stage where it can seriously 
be compared to data, this represents a revolution in the manner in which statistical physics can be applied to 
predicting observables for finite systems. We foresee that this approach will soon be developed to the stage where it 
can more conveniently, and more transparently, provide insight into the interpretation of a variety of experimental 
measurements. 



VIII. APPENDIX: DERIVATION OF RECURSIONS RELATIONS USED IN CANONICAL ENSEMBLES 



A. Recursion relation for the partition function 



The recursion relation described here was first presented by Lee and Mekjian p3| , and was first applied in a 
liquid-drop context by Chase and Mekjian Cs). 

The general relation for the partition function of non-interacting species is 

^= e n^r- ^ 

<Sfifcfl^- A> k 

where Uk is the partition function for a single particle of the species k which has size a^. For each term in the sum, 
one can factor a term uo k out of the partition function. By using the fact that Hn^aujA = 1, one may rewrite the 
partition function, 

k <Sni«i=A> K i^k 1 

k <Sn^ai— A— a k > i 

From combining this expression with Eq. (|3^), one can extract the recursion relation, 

0,A = ^2^k-jflA-a k ' (38) 



B. The recursion relation for factorial moments 



Factorial moments allow convenient calculation of the multiplicity distribution as seen in the subsequent subsection. 
Given the partition function, £7^, the moments, F^A,mi are trivial to calculate for an individual species, k. 

F k ,A, m = (rik(nk - l)(nfc - 2) ■ ■ • (n k - m + 1)) (39) 

1 uj ni 
= o~ V n fe (n fc -l)---(n fe -m+l)TT-^ r 
s I a . . nA 

<&l i es i =A> i 
_ n ^A-ma k 



J k — a • 

However, they are more difficult to obtain when they are defined in terms of n b comprised of several species with 
different masses, 

n b = ^2n k (40) 

However in this case one may proceed with the help of a recursion relation for the factorial moments. To derive the 
recursion relation, we consider the function /. 

f A ,m{b)= ]T Y[^-n b (n b -l)---(n b -m + l) (41) 

p _ fA,n(b) , . . 

fb,A,n — — o l 4 ^J 

Ua 

For the first term n b = ^2 keb n k in the sequence of m terms n b (n b — 1) • • • (n b — m + 1), each power of n k may be used 
to cancel n k in the factorial. By then factoring ujk outside the sum over configurations, one may rewrite / as 

/A,m(6) = X) w * E T\^<(n' b -l){n' h -2)---{n' b -m), (43) 

keb (En k i=A—a k i) k' 



where n' b represents the number of b- type fragments in the set k' , which differs from the previous set by the reduction 
of one fragment of type k. From the definition of /, one may rewrite Eq. ([f3|), 

fA,m{b) = / J UkfA-a h fn-\, (44) 
keb 

which leads to the recursive expression for F, 

± b,A,m — a k ,m — l ? 

(45) 

— *'a 

keb A 

This recursion relation allows one to calculate factorial moments of increasing order and for increasingly large nuclei 
given knowledge of the partition function. 



C. Obtaining the multiplicity distribution from the factorial moments 

One can express the multiplicity distribution Pb,A,n in terms of the factorial moments. Here, Pb,A,n is the probability 
of observing n fragments of type b in an event from a system of mass A. The desired expression, which we derive 
further below, has the simple form, 

Pb.A.n = E *M,m 7 ^H m ""' (46) 

m>n y 1 

where Fb t A,m = (nb(nb — 1) • • • (rib — m + 1)). Only factorial moments of order n or greater contribute to Pb,A,n since 
events with multiplicity rib < rn do not contribute to Fb^A.m- 

To prove Eq. (^6|) we rewrite the right-hand side of Eq. ( ff6| ) using the definition of factorial moments, 

E *w? — ^H ra " = E E p M^77-% 7 — 1 -^ n (-) m_n ( 4? ) 

' (to — n )!n! ' £ — to ! m — n)!n! 

m>n v 7 m>nl>m v ' v ; 

= y y Pb,Ajj- D \T — n( _ ) m ' 

— ' z — ' (£ — n — to ! TOlrt! 

where, in practice, the sums do not extend to oo due to the finite size of the system. The sum over to can now be 
eliminated by using the identity, 

E n fc \i = (i - - 4,o, (48) 

L — ' K — to !to! 

m=0 v ' 



to obtain Eq. ((46|). Although Eq. (^6|) is easy to implement numerically, it is susceptible to problems with numerical 
accuracy due to the alternating signs. Our experience is that such problems set in when the multiplicities approach 
or exceed ten. However, a recursion relation for the multiplicity distribution, which is derived in the next section, 
allows calculation of the multiplicity distribution without first calculating the moments. Such an expression does not 
have alternating signs and therefore is less susceptible to numerical problems. 



D. The recursion relation for the multiplicity distribution 

In the previous sections of the appendix, relations have been derived that give a recursion relation for the factorial 
moments, and also give the multiplicity distribution in terms of the factorial moments. In this section we derive 
a recursion relation for the multiplicity distribution, that will allow the calculation of the multiplicity distribution 
without first calculating the moments. 

By inserting the recursion relation for factorial moments, Eq. (|4^), into the formula for deriving the multiplicity 
distribution in terms of factorial moments, Eq. (|i^), one obtains, 



p, 



b,A,n — 



E? ?L W-) m ~"£ w * F M-«*.m-i%^ ( 49 ) 

^ (m — ny.nl ^ VL A 

m>n y ' keb 



f2/i ' m!n. 

feefc m>0 

By replacing -F&,A— a*,ra+m— l m the above expression with it's definition in terms of the multiplicity distribution, 
P V/ . 1 f y» V P (n + m + ro'-l)! 

fceb m>0 m'>0 ' 

^A-a fc \ - p (n_+M - 1)! 1 

— / ^fe — / -fft-A— Ofc,n+M— 1 i / l~ ) — T777 \T 

^— ' Vt A ^ n\ m\(M-m)\ 

keb A M>0 0<m<M v ; 

keb A 

where the last step utilized the identity, Eq. (flq). 

In practice, the multiplicity distributions are calculated for small A, then for successively larger A using the recursion 
relation above. However, calculation of the n = term can not be determined from the recursion relation and must 
be determined through the constraint, ^ P n = 1. 
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FIG. 1. Mass distributions are displayed (scaled by a) at three temperatures for an A — 250 system. At the high temperature, 
8.2 MeV, (dotted line) the distribution is dominated by small fragments, while below the fragmentation threshold, at 7 MeV, 
(solid line), most nucleons reside in one fragment. At 7.8 MeV (dashed line), the system seems equally divided between the 
large cluster and small clusters. The transition occurs over a remarkably narrow range of temperatures. 
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FIG. 2. The average multiplicity of IMFs, those fragments with masses from 6 to 40, are shown in the upper panel as a 
function of temperature for an A — 100 system. Including Coulomb (circles) lowers the threshold temperature for fragmentation 
compared to when Coulomb is ignored (triangles). Correlation coefficients are displayed in the lower panel. The approximations 
described in the text, which are represented by lines, describe the behavior at low excitation. A positive correlation at the 
fragmentation threshold arises from two classes of events, those with and without a large cluster, both contributing at the 
fragmentation threshold. 
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FIG. 3. The evolution of the correlation coefficient is shown for a small range of temperatures for increasing large system 
sizes. For large systems the peak narrows into a singularity, a signal of discontinuous behavior. This is remarkable given that 
it is the volume, not the pressure, that is held constant while the temperature is changed. Coulomb has been neglected in this 
calculation. 



0.3 




- 0.2 



0.1 



1$ 



l — i — i — r 



n — i — i — r 



i — i — i — r 



* 



\ a aaaa. 



4/ 



A <^ 



0.0 ^^o^oWBd^^DaDaDaDaD'fe 



o 



10 



20 



30 



n 



IMF 



FIG. 4. Multiplicity distributions are shown at three temperatures for the A — 250 system. Below the fragmentation 
threshold at a temperature of 7 MeV (squares) , the distribution is peaked at a low multiplicity and is well described by a negative 
binomial distribution shown as a dashed line. At a high temperature of 9 MeV (circles), the distribution is sub-Poissonian 
and well described by the binomial distribution (dashed line). At 8 MeV (triangles), the multiplicity distribution is strongly 
super-Poissonian and is not even well described by a negative binomial. 
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FIG. 5. IMF multiplicity distributions are displayed in the upper panel, while correlation coefficients are displayed in the 
lower panel. Percolation calculations, based on a spherical lattice of 203 sites are represented by triangles while the canonical 
ensemble, with A — 200 and no Coulomb effects, is represented by circles. The percolation model does not exhibit a first order 
phase transition, hence it's correlation coefficient does not show a spike at the critical temperature. Other notable differences 
are that the multiplicity in the statistical calculation reaches higher values than in the percolation model while the correlation 
coefficient in the percolation calculation reaches larger values at low temperatures. 



